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NONLINEAR ANALYSIS AND PERFORMANCE EVALUATION OF THE 
ANNULAR SUSPENSION AND POINTING SYSTEM (ASPS) 

By 

Suresh M. Joshi* 


SUMMARY 


The Annular Suspension and Pointing System (ASPS) can provide highly 
accurate fine pointing for a variety of solar-, stellar-, and Earth-viewing 
missions. In this report, a detailed nonlinear mathematical model is 
developed for the ASPS/Space Shuttle system. The equations are augmented 
with models of components such as magnetic actuators and gimbal torquers. 
Control systems and payload attitude state estimators are designed in 
order to obtain satisfactory pointing performance, and statistical pointing 
performance is predicted in the presence of measurement noise and 
disturbances. 

INTRODUCTION 


The Annular Suspension and Pointing System (ASPS) holds the promise 
of providing very high payload pointing accuracy and stability during 
space shuttle orbital missions. The basic concept of ASPS was suggested 
in reference 1, and was developed in reference. 2, along with preliminary 
results on statistical pointing performance based on a linearized mathe- 
matical model. A more elaborate linear analysis was carried out in 
references 3 and 4 in which a detailed model of the ASPS/Space Shuttle 
system was developed and used. Payload attitude measurement noise and 
state estimators were also included in references. 3 and 4. The linear 
analyses used linear mathematical models for stochastic crew motion 
disturbances. 

* Research Assistant Professor, School of Engineering, Old Dominion Univer- 
sity, Norfolk, Virginia 23508. 



The physical system, however, contains a number of nonlinearities. 

The nonlinearities include inherent trigonometric nonlinearities, magnetic 
actuator nonlinearities, gimbal friction and torquer nonlinearities, etc. 

It was therefore necessary to include these nonlinearities in order to 
more accurately predict the performance of ASPS. Detailed information 
regarding ASPS components (e.g., magnetic actuator models) also became 
available from the Sperry Flight Systems, under contract NAS1-14214, 
during the course of this work, and was used in the nonlinear mathematical 
model . 

The ASPS (fig. 1) consists of an elevation gimbal and a lateral gimbal 
for "coarse" pointing. Magnetic suspension and fine pointing actuators 
are provided in the vernier assembly and generate fine pointing control 
torques. The payload or vernier assembly (VA) , which consists of an 
annular rim (made of magnetic material) attached to a payload mounting 
plate (PMP), is levitated by the magnetic actuators of the VA. The payload 
instrument is mounted on the PMP. A noncontacting optical data coupler 
is used for transmission of data between VA and the Space Shuttle pallet. 

A battery pack mounted on the PMP supplies the payload power and eliminates 
the necessity of any cables. As a result of design evolution, the present 
version of ASPS, which is shown in figure 1, is somewhat different from 
the ASPS used in references 3 and 4 (although mathematical descriptions 
do not significantly differ) . The functions of various components remain 
the same, and are described in references 2, 3, and 4. The magnetic 
actuator configuration is also different from that used in references 2, 

3, and 4. There are three axial actuator stations (A, B, C) spaced 
120° apart and two radial actuator stations (U,V) spaced 90° apart, as 
shown in figure 2. Roll motor station W consists of two segments, one 
inside and one outside the rim. The system considered in this study has 
been sized for payloads up to 600 kg with Z-axis center of mass offsets 
up to 1.5 m. The magnetic suspension actuators use bias current lineariza- 
tion to remove the current squared nonlinearity. To compensate for the 
inverse-gap- squared relationship, a signal proportional to the gap is 
used to multiply coil currents. A segmented two-phase solid iron rotor 
AC induction motor controls the roll rotational servo. Proximity sensors 
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are provided at each axial and radial actuator station, as well as at the 
roll motor station to compensate for radial gap. The "coarse" gimbal 
torquers are permanent magnet, two-phase brushless DC type. 

The objectives of this study are 

• Development of complete nonlinear equations of motion for 
ASPS/Space Shuttle system. 

• Inclusion of all available component models (e.g., magnetic 
actuators, "coarse" gimbal torquers, etc.). 

• Design of control systems for the vernier assembly and the 
"coarse" gimbals. 

• Design of payload attitude state estimators. 

• Prediction of statistical pointing performance via digital 
computer simulation. 

The analysis presented here is based on rigid-body models, and thus 
does not contain representations of structural modes. 

SYMBOLS 


A, B, C 
A, A, A , 


A. . 
il 


A- 

1 


B. 

1 

C C 

V 

c 

g 

C (x) 
p 

D. . 
il 

d 

E 


Axial actuator stations 
Coefficient matrices 

Point on lateral gimbal representing actuator station 
(j = A, B, C, U, V) 

Point on vernier assembly corresponding to A^ 
Coefficient matrices used in controller design 
"Coarse" gimbal control matrix 
Cross product matrix of x 

Transformation matrix from i- to j -coordinate system 
Five-dimensional vector appearing in eq. (67) 

Body rate to Euler rate transformation matrix 



e . 
1 


F 

F . 
ai 

F bias 

F 


ca 

F 
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Unit vector in i direction (i = 1, 2, 3) 

Magnetic actuator force vector (total) 

Magnetic actuator axial forces 
Cable bias force 

Total cable force 

Cable spring force 

Force vector generated by actuator station j 

Force exerted by lateral gimbal on elevation gimbal 
contact surface 

Reactive force acting on coarse assembly 
Radial disturbance force due to roll motor 
Magnetic actuator radial forces 
Force acting on shuttle 
Tangential force due to roll motor 
Force vector defined in eq. (78) 

Cogging frequency 

Ripple frequency 

Inertia matrix of shuttle 

Inertia matrix of elevation gimbal 

Inertia matrix of lateral gimbal 

Inertia matrix of vernier (payload) assembly 

Performance function 

Control gain 

• • 

Vector functions of u,<l>,d, <p, 9, a, a, u 

o o v , v 

Proportional gain for gimbal control system 
External moment acting on shuttle 

Magnetic actuator torque acting on vernier assembly 



m 


m 0 


mi 


m2 

m 

v 

0 


0. 

1 


0 


P 

P 


cmi 


cm 2 


cma 


R 


cmm 


R £ei 


r 

cmi 


m 0 + mi + m 2 

Mass of shuttle 

Mass of elevation gimbal 


Mass of lateral gimbal 
Mass of vernier assembly 

Origin of 2"-frame; also the point fixed to lateral 
gimbal which is the nominal position of payload 
mounting plate center 

Origin of i -frame 

Center of payload mounting plate 

Center of mass of shuttle 

Combined center of mass of shuttle and "coarse" gimbals 


Center of elevation gimbal axis 
Center of lateral gimbal axis 
Percent ripple 

Position of elevation gimbal c.m. relative to P in 
0-system 

Position of lateral gimbal . c .m . relative to P^ in 
0-system 

Position of lateral gimbal c.m. relative to P g 
in 1-system 


mi 



+ ro 2 R 


cm 2 


Position of a point on the mutual contact surface 
between lateral and elevation gimbals relative to P g 


Radius of payload mounting plate used in magnetic 
actuator torque expressions 


Position of combined c.nr. of, shuttle and two gimbals 

relative to ; 0. 

1 


Position of elevation gimbal c.m. relative to P g 
in 1 ' -system - 
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r 

cm 2 


r 

cmv 

r d 

r 

e 


r . 
J 


r 


% 



r 

P 


r l> r 2 


rj, r 


(f)VC 

T . 
bias 


ca 


cog 


cogo 


com 


fo 


fil 


T. 

1 


t 

2 


Position of lateral gimbal c.m. relative to P in 
2 '-system 1 

Position of vernier assembly c.m. in v-system 

Position of P relative to P in 1' -system 

Position of payload mounting plate center relative to 
0 in 2"-system 

Position of actuator station j in 2"-system 

Position of 0 relative to lateral c.m. in 
2-system 

Position of a point on the mutual contact surface between 
lateral and elevation gimbals relative to . P 

Position of P relative to 0. 

i 

Position of vernier assembly c.m. relative to 

0 . 

l 

Positions of c.m. of elevation and lateral gimbals 

relative to 0. 

i 

Positions of P and P„ relative to 0. 

e Z i 

Rate command for . d> 

Y v 

Bias cable torque 

Cable torque acting on vernier assembly 

Cable torque acting on lateral gimbal 

Cogging torque 

Peak cogging torque 

Commanded torque 

Disturbance torque 

Gimbal friction torque 

Peak or Coulomb friction torque 

Estimator update interval 

Noise-generation time interval 
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T 

RCJ 

T 

r 

T • , 
ripple 

T l 

T l 

^com 

t 2 

T2 com 

T i(v) 

T 

V 

T 3 

U, V 

(X, Y, Z) 


z 

a . 
i 

a 

o 

a 


v 


a 


vr 


3 

3 

c 

e i 

r l> r 2 
Y 

AA, AB, etc. 


6. 

J 


Reaction control jet torque 

Reaction torque acting on lateral gimbal 

Ripple torque 

Elevation gimbal torquer torque 
Elevation gimbal commanded torque 
Lateral gimbal torquer torque - 
Lateral gimbal commanded torque 

Transformation matrix for v radians rotation about 
axis i 

Magnetic torque acting on vernier assembly 
T 3 (45°) 

Radial actuator stations 

Denotes coordinate frames; subscripts are described 
in the text 

Integrator state variable 
Coefficient matrices 
Attitude vector of shuttle 
Attitude vector of vernier assembly 

Attitude vector of vernier assembly relative to lateral 
gimbal 

Incremental gimbal angle vector 

Commanded gimbal angle vector 

Vector functions of state variables 

Matrices used in payload following controller design 

Gimbal torquer angle 

Linearized center errors and gimbal angles [defined 
in eq. (67)] 

Centering error vector at actuator station j 
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Desired decay characteristic 
Damping ratio 

Position of actuator station j relative to 0 

Position of point on rim corresponding to A. 
relative to 0 


Position of combined c.m. of shuttle and two gimbals 
relative to P 


cmv 


cmi 


cm2 


Pl> P2 


Position of vernier assembly c.m. relative to 0^ 

Position of elevation gimbal c.m. relative to P g 

Position of lateral gimbal c.m. relative to P^ 

Position of P. relative to P 
£ e 

Position of 0 relative to 0 

v ' 

Position of 0 relative to lateral gimbal c.m. 

Positions of elevation and lateral gimbal c.m. relative 
to P 


! 1 

Pi) P 2 


se 


Positions of P and P„ relative to P 
e £ 

Torque exerted by shuttle on elevation gimbal 


(<f> , 

0) 

Elevation and lateral gimbal angles 


0 , ip ) 
o’ r o 

Components of a Q 

<V 

0 , Ip ) 
V r v 

Components of 

^vc 5 

0 , lj> 

VC VC 

Commanded values of <(> , 0 y , 

(<t> 

^ Y vr 

, 0 , <J> ) 

vr vr 

Components of a 
r vr 

n 

VT 


Angular velocity vector of vernier assembly relative 

V J. 


to lateral gimbal 



Angular velocity of elevation gimbal relative to 
shuttle 

9-2 


Angular velocity of lateral gimbal relative to 
elevation gimbal 

w 


Natural frequency 
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d) 0 , mi) u >2 Angular velocities of shuttle, lateral, and elevation 

gimbals relative to i-frame 

a) Angular velocity of vernier assembly relative to 

v i-frame 

Angular velocity of 1' -frame 

Angular velocity of 2' -system 

Position of actuator station j in v-system 

An overhead bar denotes a vector quantity and corresponding operations 
indicate vector operations. (Overhead bars are not used to denote vectors 
in the Euclidian sense.) Superscripts "T” and "-1" denote transpose 
and inverse, respectively. One and two overhead dots denote the first 
and second time derivatives, and an overhead ,,A " denotes the estimated 
value of a variable. 

MATHEMATICAL MODEL DEVELOPMENT 
Coordinate System 

Referring to figure 3 (xppp is an inertial coordinate system centered 
at CL. (XqYqZo) is a shuttle-fixed system centered at P, the center of 
mass of the shuttle. (XJYJZJ) is a system fixed to the elevation gimbal, 
with origin at P , the center of the elevation gimbal rotation axis. 

It is obtained from (X 0 Y 0 Zo) by one positive rotation <j> (elevation gimbal 
angle) about the X 0 axis. Coordinate system (XiY^)) is also fixed to 
the elevation gimbal and is parallel to the (XjYJZJ) system. It is cen- 
tered at the elevation gimbal center of mass. System (X^Y^Z^) is obtained 
by one positive rotation 0 (lateral gimbal angle) about the Yj axis, 
and is fixed to the lateral gimbal. Its origin is the center, P^, of 
the gimbal rotation axis. System (X 2 Y 2 Z 2 ) , which is also fixed to the 
lateral gimbal, is parallel to (Xppp and has the lateral gimbal center 
of mass as its origin. (X^'Y^'Zp is a system fixed to the lateral gimbal 
which is parallel to (X 2 Y 2 Z 2 ) • Its origin lies at 0, the center of the 
lateral gimbal plate, which is also the nominal position of the center (Cp 
of the payload mounting plate (PMP) . (X y Y y Z ) is a coordinate system fixed 


“1 

“2b 

jv 
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to the vernier assembly (VA) , and its origin is the center, 0 , of the 
PMP . 

* 

"Coarse" Assembly Equations 

The "coarse" assembly consists of space shuttle, elevation gimbal, 
and lateral gimbal. Equations of motion for the "coarse" assembly can 
be derived by applying the method given in reference 5 (pp. 374-377) for 
a system of rigid bodies. "Coarse" assembly translation is given by: 


m r 

c 


F 


r 


+ F 

s 


( 1 ) 


where m = mg + m^ + m 2 , and r £ is the position of the combined c.m. of 
the space shuttle, elevation gimbal, and lateral gimbal relative to Ch. 

F^ is the total reactive force acting on the "coarse" assembly, and F s 
is the sum of other forces acting on the shuttle. F^ consists of magnetic 
actuator force and cable force. The latter is present only when a cable 
is used for data transmission as a backup for the optical data coupler. 
Since 



and 

p c = m *- mi Pl + m2 P2 ^ 
equation (1) becomes 

t 

m r p + m l Pl + m2 P2 = F r + F s 

(see fig. 3 for explanation). 

"Coarse" assembly rotational equations are given by 

H = M - p xmr 
p p c p 


( 2 ) 


(3) 


(4) 


( 5 ) 
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where is the angular momentum with P as the reference point or 

• • 

Hp = 1q ujq + Ij o) i + I 2 ( 1)2 + Pi x mj pj + P 2 x m 2 P 2 (6) 

where Ii and I 2 are the inertia matrices of the elevation and lateral 
gimbals relative to the 1- and 2-system, as observed by an observer 
rotating with the 0-frame. That is 

*1 = Do 1 I 1 B D 01 > l 2 = d o2 d o2 ( 7 ) 

where I 1d and l 2 n are inertia matrices of elevation and lateral gimbals 
relative to the 1-frame and 2-frame, respectively. is the torque 

exerted on the "coarse" assembly. 



T 

0)1 - Wo + D 0 i 51i - wq + f^i 

(8) 


T T 

W 2 = 0)o + + Dq2 ^2 - ^0 + ^1 + Dq 1 f^2 

(9) 

Using 

the notation 



H. = T. w. 
1 11 



• # • • 

H = Ho + Hi + H 2 + Pi x mi Pi + P 2 x m 2 P 2 = M - P x m 
P P c 

t (10) 


• • 

Hi = (H i ) r + w 0 x H i , i = 0,1,2 

(ID 

where 

CH.) is the derivative of H. relative to the 0-frame, 

a/r 1 

expressed 

along 

0-coordinate system, 



• 

(H.) =1. w. +I.o). , i = 0,1,2 

1 r 1 1 1 1 

(12) 


• • 

p i = (P i ) r + 2w 0 x (p i ) r + w 0 x (w 0 x p i ) + w 0 x p i5 i = 

1,2 (13) 
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But 


Pi = Pi + P 


cmi * 


p 2 = Pi + P d + P 


cm 2 


(14) 


and 


(P l) r = (Pl) r = o 


(pi) = pi + dL r A R 
1 r 1 Ui cmi = cmi 


(15) 


^ p d^r r c 


— i T T 

(p?) = Pi + Dm r, + Dn 2 r A R 

uid uz cm2 = cm2 


(16) 


Components of vector derivatives pj, p 2 can therefore be evaluated as in 
equation (13) along the 0-axes to obtain 


Pl - dj (0 q + d 2 (J> + $1 


(17) 


p 2 - ot 3 ojq + <j> + 0 + f3 2 


(18) 


where matrices ou, 8^ are functions of co 0 , <f>» S, and 0. 

"Coarse" assembly rotational equations of motion are then obtained by 
substituting for m r^ from equation (4), and for pi and p 2 from 
equations (17) and (18), into equation (10): 


Aii (On + A 12 <|) + A13 0 + L -M - — C (R ) (Dq 2 F + F ) 
u r i;S mo pm p v cmm r s 


(19) 


where matrices A. . and L are functions of eon, • <f>, 0, 6, 0, and 

ij u>o u 


R = mi R + m 2 R 
cmm 1 cmi cm 2 


( 20 ) 
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To obtain complete equations of motion for the "coarse" assembly, 
equations for elevation and lateral gimbals must also be derived. In 
these derivations, forces acting at the contact surfaces between shuttle 
and elevation gimbal, and between elevation and lateral gimbals, must be 
considered. Forces exerted by the elevation gimbal on the lateral gimbal 
at their mutual contact surface are assumed to be acting radially, i.e., 
they produce no torque about the lateral gimbal axis of rotation. Similarly, 
forces exerted by the shuttle on the elevation gimbal at their mutual 

t 

contact surface are assumed to produce no torque about the X^-axis. 

Equal and opposite reactive forces act on corresponding bodies at contact 
surfaces. Since these internal forces must be eliminated from the equations 
of motion, it is convenient to derive lateral gimbal equations first. Let 
F^ei be a force exerted by the lateral gimbal on the elevation gimbal 
at point (with respect to 2' -system) located on the contact surface 

(i = 1,2,3...). Force -F^ ^ is exerted by the elevation gimbal on the 
lateral gimbal. 

Lateral gimbal translation equation is 

■2 - F T - £ F-*! (21) 

But 

r 2 = r p + p" 2 (22) 


Substituting for r 2 into equation (21), and for r p from equation (4), 


Ef 


fcei 


T&2 \ ■ mi m 2 •• 

= (1 F + pj 

m / r . m 1 



1 P2 



(23) 
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It is convenient to derive lateral gimbal rotation equations in the 
2' -system. Let I 2 g be the inertia matrix of the lateral gimbal relative 
to the 2' -system. Then 

H 2 = I 2 g w 2g C 24 ) 

where 

u 2g = D 02 “0 + D 12 ^1 + ^2 

(fi 2 ) r + -2g x (H 2 ) r = T - 7 £e . x F £e . + f 2 

where is the reaction torque (consisting of magnetic torque and cable 

torque), and T 2 is the torquer torque. r 2 is the position vector of 

P. relative to 0. . 

£ i 


(25) 


P X mo To 

cm 2 z z 


(26) 


r 2 =r +p 2 -p 
z P cm 2 


(27) 


After evaluating r 2 . in a way similar to equation (13), and simplifying 


a 31 w 0 + a 32 


+ a 33 8 + l' 


- T r * t 2 - (E 


r . x F. . 
£ei £ei 


Hl2 


C (r ) (F + D 02 F ) 
m p cm 2 r s 


(28) 


Equation (28) is obtained along the 2’-system. Matrices A. . and L Q 

are functions of wq, <f>, 8, <P, and 8. Equation (28) gives the value 

of r . x F„ . (along the 2 ' -coordinate system). Since forces F„ 

£ei £ei 6 , _ f . £ei 

do not produce torque about the Y 2 axis, the Y 2 component of 
^ r £ei x ^£ei '*' s zer ° - Thus the lateral gimbal equation is 
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( 29 ) 


A31 “0 + A3 2. 4 > + A3 3 6 + L 0 = ?2 ■ (T r + T 2 ) 


m 2 T> 

— e C (r ) (F + D 02 F ) 
m p cm2 r s 


where 


a T . 1 t T . 1 

A 3 j - e 2 A 3 j , Lg ~ e 2 L 


0 


Elevation gimbal rotation equations are conveniently written in the 
1 ' -system. Let 

H i = ii B “i B 

I ; . ' , 

where is the inertia matrix of the elevation gimbal relative to the 

1 ' -system. 

“lg = ° 01 w ° + ^ 


CH'D r > “i B X Hi * R tei * F tei * T, * T se - p cmi X m, r| - T 2 


where R . is the point of application of force F . exerted by the 
lateral gimbal at the contact surface, and x se is the torque due to 
forces exerted by the shuttle at the elevation gimbal-shuttle contact 
surface. Tj is the torque exerted by the torquer, and -T2 is the 
reaction torque exerted by the lateral gimbal torquer 


R £ei P d + r £ei 


(Hi) r + «! b x H X = p d x.F 4e . + £ r £ei X + T 1 


- p x mi ri + x - To 
cmi 11 se 


( 30 ) 

( 31 ) 

( 32 ) 


( 33 ) 

( 34 ) 
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x se has no component along the Xi axis 


r l = r p + p l 


cmi 


( 35 ) 


The elevation gimbal equation can be obtained by substituting from equation 
(28) for r £ e i x F £ e i’ anc * ^ rom equation (23) for ^ F a ^ ter evaluating 
P cm x mi r{ as in the case of the lateral gimbal, and taking only the 
first components 


r 


A 2 i -u) 0 + A 22 <t> + A 2 3 0 + = ei 


Ti + Dl 2 T 


/mr 


- C ( — R + — r ) 
p \ m cma m cmi/ 


( D 


12 F r + D 01 F s 


) 


+ C (r.) D 2 F 

p>- QJ iz r 


] 


(36) 


where 


R 


cma 


r d + ° 12 r 


cm 2 


(37) 


Vector coefficients A., and L are functions of wq, 4»., 6, 4>, and 
0. Thus equations (19), (29), and (36) define complete rotational motion 
of the "coarse" assembly. 


Vernier Assembly Equations 


Let r^ be the position of vernier assembly (VA) center of mass relative 

to 0. . Then 
i 


m 


r = 



(38) 
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But 


r 

v 


r 

P 


+ P2 + 


+ P 

cmv 


(39) 


p can be evaluated as in equation (13). It is convenient to evaluate it 

X/ ■■ 

along the 2 ' -coordinates . P cmv can conveniently evaluated along the 
v-system: 


p = (i) x p 
cmv v cmv 


0) X 
V 


(a) x p ) 
v v cmv 


(40) 


p , the acceleration of PMP center relative to 2"-frame, can be evaluated 
e’ 

along 2"-coordinates as in equation (13). Substitution and simplification 
yield the following VA translation equation along the 2"-system: 


r + 
e 


H 11 w 0 + H I2 <t> + H 13 9 + H 14 


0) 


+ 


L 

v 



- D 02 F 
m s 


(41) 


where and L v are functions of mo, 

and r .. 
e 

Defining 





“0* < f'» 9 « v » 


H = I to 

V V V 


(42) 


where I is the VA inertia matrix with respect to the v-system 


I to + o) x (I a) ) = T + T - p xm r 
vv v v v v c cmv v v 


(43) 


where T v and T c are magnetic actuator and cable torques in v-system. 
Substituting for m v r v from equation (38) and expressing along 2"- 
coordinates results in. 
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w = I -1 -C (to ) I to + T - C(r ) 
v v p v v v v p cmv 


x (F - Do F ) + T 
v z v ca ; c 


] 


( 44 ) 


Combined Equations of Motion 
The "coarse" assembly equations have the form: 


<*>0 

<f> 

0 


= y 


(45) 


where A is a 5 x 5 matrix and y is a 5 x 1 vector of functions of 

• • 

^0* 6, <(>, and 0, and the forcing functions. These equations 

can be solved for (wq, <j>, 0) by inverting A: 


•T •• T -i 
(o) 0j <f>, 0) = A 1 y 


(46) 


Other equations are 


«0 = Eq 1 w o 


(47) 


dt + = * 


(48) 


dt 


0 = 0 


(49) 


dt r e r e 


a = E" 1 ^ 
v v v 


(50) 

(51) 


where Eg and E y are Euler to body rate transformation matrices for 
shuttle and VA. 
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Equations (46) to (51) , along with (41) and (44) represent complete 
rigid-body equations of motion for the ASPS/space shuttle. 

Magnetic Actuator Centering Errors 

Let denote the position of actuator station j relative to 0, 

the lateral gimbal plate center, and Pg^ denote the position of corre- 
sponding point Bj in the vernier assembly, relative to 0^. Then the 
magnetic actuator centering errors are given by 


6 = -p . + p + p„ . 

1 p Aj p e p Bj 


In the 2"-system, this equation reduces to: 


(6.) = r + (D 0 p D-o d! - I)r . 

v j'r e i u iv j 


where r^ = position of actuator station j in 2"-system, and 


D iv = T 2 (0 v ) Tl (* y ) 


For axial actuator stations. 


-/3/2 


r A = r 1 , r B = r -1/2 , r c = r I -1/2 


For radial stations. 


- 1//2 


r = r 1/J2 , r = r 

u v 


1//2 

1//2 
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For roll motor, 

0 
-1 
0 

At axial actuators, only axial centering errors are measured, which are 
given by the third components of 6^, 6^, 6^. Radial centering errors 

at stations U and V are obtained by transforming 6^, 6^ by 45 degrees 

rotation about the roll axis and then taking appropriate components. Axial 
and radial centering error at station W are given by third and second 
components of 6^. 

It is also necessary to compute the time derivatives of centering 
errors for use in the magnetic actuator models. They are computed along 
the 2"~system using the method of equation (13). 

Magnetic Actuator Models 

The radial and axial magnetic actuators, which provide rim suspension 
forces as well as fine pointing torques, are inherently nonlinear devices. 
The nonlinearities are present because of the fact that the magnetic force 
is proportional to current squared, and is inversely proportional to the 
square of the gap. A bias current technique has been used to remove the 
current-squared nonlinearity, and a signal proportional to the measured 
gap is used to multiply each coil current in order to remove the inverse 
gap squared nonlinearity. 

A mathematical model of such a magnetic actuator (axial or radial) was 
developed by Sperry Flight Systems under contract NASI- 14214. Each actuator 
model is of order 6 (3 for each of the 2 electromagnets) . The model also 
incorporates limiters and losses due to eddy currents. 

Roll Motor Disturbance Model 

The roll motor selected for generating tangential force at station W 
is an ac induction motor with constant fixed field excitation and a control 
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field excitation which is proportional to the tangential force command. It 
consists of two curved stator segments, one on the inside and the other on 
the outside of the rim. The rim is nominally equidistant from both segments. 
Each segment has a fixed and a control field winding. Since both segments 
exert attractive forces on the rim which are inversely proportional to the 
gap squared, each stator winding is made to carry a current which varies 
linearly with gap in order to achieve radial force balance. A small magni- 
tude anomaly radial force is still present. The roll motor is excited by 
a 100-Hz carrier--therefore, it is not practical to include a complete roll 
motor model in view of the enormous computer time requirements. However, 
a "worst" case situation can be simulated by including a 200-Hz (twice the 
carrier frequency) sinusoidal radial disturbance force at the roll motor 
station. A 2- or 3-second time history should then yield adequate results 
indicating effect of the roll motor. 


Cable Forces and Torques 

During normal operation, a noncontacting optical data coupler is used 
for transmission of information between the payload and space shuttle 
pallet. However, a transmission cable may be used as a backup system for 
the data coupler. The cable is assumed to be routed through the lateral 
"coarse" gimbal and the center of payload mounting plate, and is assumed 
to cause a constant axial force, a constant radial force, a constant 
torque about the roll axis, as well as axial and radial spring forces 
(proportional to relative displacements) and a torsion torque (proportional 
to relative roll angle). This simple representation is used to keep the 
number of state. variables low.. In 2"-coordinates , the spring forces acting 
on the lateral gimbal are given by 


F . = K . r . , i = x,y,z 
csi si ei J 


(58) 


The total cable force acting on the lateral gimbal (in 2"-coordinates) is 
given by 


F - F + F, . 
ca cs bias 
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Relative Euler angles between the payload and the lateral gimbal are 
computed from: 


a 


vr 


= E _1 a 
vr vr 


( 59 ) 


= D 0) 009 n 

vr v 2 v 

T 

where a ^ = (<j> , 8 vr , ^ vr ) is the attitude of the payload (Euler angles) 
relative to the lateral gimbal, and E“^ is the body rate to Euler rate 
transformation matrix. The cable torque acting on the lateral gimbal is 
given by: 



0 


0 

T 

0 

+ 

0 

ca 





_ T bias 


_ k ts V 


( 60 ) 


The cable torque acting on the VA is given by: 


T = D 2 T 
c ^v ca 


( 61 ) 


where T £ is expressed in the v-system, and D 2 ^ denotes the transformation 
matrix from the 2- to v-system. 

Gimbal Torquer Models 

A permanent magnet brushless dc torque motor is used for driving each 
"coarse" gimbal. It is assumed that the current loop has been compensated 
to obtain a high bandwidth, so that 


T=T -T^ + T + T . . 

out com f cog ripple 


( 62 ) 


where T ., T , T- are the actual, commanded, and friction torques, 

out’ com’ f ’ ’ n > 

and T co g and ^ripple are c0 §8i n § anc * r ippl e torques. Friction 

torque T^ is generated using a Dahl model (ref. 6) 
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dY 


fo 


Sgn(y) 


where 

Y = girabal angle (rad) 

= rest stiffness or slope (NM/rad) 

= peak or Coulomb friction torque (NM) 
i = the friction law (shape of the hysteresis curve) . 
The friction torque is given by 


dT. dT- , 
f _ f dY 

dt dY dt 


Ripple torque is given by: 


T 

! ripple 


* Pr T com (l Sin ( f r ^l ‘ *) 


where 

p r = percent of ripple 

f r = ripple frequency = no. of phases x no. of pole pairs. 
Cogging torque is given by: 


T = T 
cog cogo 


sin 


(f c y) 


where 


T cogo = peak c °ss in s tor q ue 
f = cogging frequency = no. of slots. 


( 63 ) 


( 64 ) 


( 65 ) 


( 66 ) 
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CONTROL SYSTEMS DESIGN 

The two control systems to be designed are the "coarse" gimbal control 
system and the magnetic suspension and fine pointing control system. The 
former may use magnetic actuator centering errors, or fine pointing error 
and rate estimates, or gimbal angle information, depending on the mode of 
operation. The magnetic suspension control system uses estimates of fine 
pointing errors and rates, and magnetic actuator centering errors measured 
by proximity sensors. 

"Coarse" Gimbal Control System 

For stellar or solar pointing missions, the basic modes of operation 
of the "coarse" gimbal control system are as follows: 

1. Fine-pointing mode, fixed gimbals. In this mode, the "coarse" 
gimbals are required to be in fixed positions relative to the Shuttle. 

2. Fine-pointing with payload following gimbals. In this mode, 
the gimbals are required to continuously follow the payload in such a 
manner as to keep the rim properly centered in the magnetic actuators. 

3. "Coarse" gimbal backup mode. In this mode, the vernier assembly 
is mechanically latched to the lateral gimbal, and gimbal command torques 
are generated using pay-load attitude and rate estimates. 

4. Gimbal slew mode. In this mode, appropriate slew rates are used 
to drive the gimbals, with the vernier assembly either latched or centered 
in the magnetic actuators . 

The generation of commanded gimbal torquer torques is quite straight- 
forward for mode 1 and mode 3. For the slew mode, the problem lies in 
computing the desired slew rates. The payload following gimbal control 
system of mode 2 can be obtained by generating gimbal angle commands 
which will minimize a norm of the axial and radial centering errors. 

This is done by linearization about the target pointing angles. The 
linearized expressions for axial and radial centering errors are 


24 



r i 


" T T _ 

< 

< 


e 3 r l r A e 3 r 2 r A 



T „ T „ 

AB 

z 


e 3 r l r B e 3 r 2 r g 



T „ T „ 

AC 

z 


e 3 r i r c e 3 r 2 r c 



T i T » „ 

AU 

r 


e 2~ T 3 r l r u e 2 t 3 r 2 r^ 



T ' T » „ 

i — 

> 

< 


e l T 3 r l r v e l T 3 r 2 r v 



+ d 


( 67 ) 


where AA , AB^, AC denote incremental axial centering errors at actuator 
stations A, B, and C, and AU^, AV^ denote incremental radial cen- 
tering errors at actuator stations U and V. A0, A<J> represent incre- 
mental lateral and elevation gimbal angles, and d is a five-dimensional 
vector. 
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r 2 = 


D 12 
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(69) 


The transformation matrices and their partial derivatives are evaluated at 
nominal values and target pointing angles. The 5x1 vector d is not a 
function of A0 and A<j>. Equation (67) can be symbolically written as: 


A = A 8 + d 
g 


(70) 


where 



is a 5 x 2 matrix, 


and B = (A8, A<p ) 


T 


Minimizing the quantity 



(71) 
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with respect to 3, the command angle is obtained as: 


3 = -CA T A ) -1 A T (A - A 3) 

c g g g^ g 


= C A + 3 
g 


Therefore, the error signal going into the control system is 


B - 


3 = C A 
g 


where 


C = - (A T A ) “ 1 A T 
g g g g 


( 72 ) 


(73) 


(74) 


The gimbal control matrix is dependent on nominal values of 

pointing angles and requires periodic updating at large intervals of 
time. Payload following gimbal control can thus be achieved by feedback 
of five proximity sensor outputs. The desired gimbal torquer torques can 
be written as : 





(75) 


is obtained from 


where 
C 

be derived from the proximity sensor outputs using second order networks. 


K , L are rate and proportional gains, and C 
gi g^ § 

by interchanging the two rows. The rate signals in equation (75) can 


Magnetic Actuator Control System 

The magnetic actuator system consists of three axial and two radial 
actuator stations and one roll motor station. Figure 2 shows the magnetic 
actuator configuration consisting of three axial actuators spaced 120 
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degrees apart, two radial actuators spaced 90 degrees apart, and roll motor 

segments producing a tangential force at one point. The axial actuator 

stations are labeled A, B, C, and the radial and roll actuator stations 

are labeled U, V, W. They produce forces F ., F ., F . The total 

3.1 n T 

forces in X 2 , Y 2 , Z 2 directions produced by the actuator stations on 

the vernier assembly are 



where 


- — — 0 0 0 1 

/2 /2 



— — 0 0 0 0 

/2 /2 


and 


0 0 111 


f = (F , F , F , F , F , F_) 
rT r 2 a l a 2 a 3 ? 


(76) 


(77) 


(78) 


Nominal torques (excluding anomaly torques) produced by these forces on the 
vernier assembly about the v-system: 



M = C x F + t; x F n + x F„ + x F„ 
v Av A Bv B Cv C Uv U 


+ ? W X F V + %v X F W 


(79) 
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where ■£ is the position of station j in the v-system, and denotes 

the force vector at the jth actuator station, and so forth. Considering 
linearized design, the desired response can be obtained by making 


M 


C (r ) F 
p cimr 


-I A 
vo av 


(80) 


where 


A = (A, , A , A. ) 
av ■ v 8v ipv 


A, = 2pw(4> - r. ) + u) 2 (<}> - 4 ) + K t z, 

<j>v r ' r v tj>vc v vc 14 >v 4>v 


z, = * - d> 

())V . T v Y vc 

Cd> , r ± , K T , , z, are attitude and rate command, integral gain, and 
vc <pvc 1 <pv (pV 

integrator output.) Or 


C T f - C (r )C f = -I A 
T p v cimr F vo av 


(81) 


where 


C T = r 


0 0 1 


0 0 0 Af- 
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-7 0 


4 • 


0 0 0 0 0 1 


(82) 


Therefore, desired pointing control can be obtained by making 

(c T - C (r )Cjf = -I A 
\ T p v cmvo F/ vo c 


av 


( 83 ) 


The desired VA translation response can be obtained by making 
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( 84 ) 


(_L + I\ C f= _ X 

\m v my F re 

where X is defined similarly to X 

re J av 

Equations (83) and (84) form a system of six linear equations in six 
unknown components of f. Therefore, the unique solution can be readily 
obtained by inverting the coefficient matrix. 

This design offers complete decoupling of VA rotation and translation, 
if the VA inertia and c.m. location are accurately known. The commanded 
magnetic actuator forces are computed in this manner, using proximity 
sensor outputs, and payload attitude and rate estimates. Magnetic 
actuator centering error rates are generated from proximity sensor outputs 
using second order networks. These force commands are fed into the magnetic 
actuators. In the simulation, the actual output forces are obtained using 
actuator mathematical models mentioned earlier. Magnetic actuator forces 
so obtained are in the 2"-system and act at points on the vernier assembly 
rim. Since the rim is continuous, only the relative roll angle tp re ^ of 
the payload assembly, relative to the lateral gimbal, needs to be considered 
for obtaining points of application of magnetic forces in v-system. The 
location of actuator station j in v-system is given by 



T 3 ( -^rel ) 


r . 
J 


- (85) 


where r^ is the position of station j in the 2"-system, which is a 
fixed vector. The three-dimensional actuator force at station j in the 
2"-system is 


F . 


l F d 


For example, for station A, F = F = 0, and F . = F , which is the 

XA y A Z A cl j 

third component of actual magnetic force vector f [corresponding to f 


29 



of eq. (78)]. For roll motor station W, the radial disturbance force (F d ) 
can be included in this formulation. The total torque on the payload 
assembly due to magnetic actuators (in v-system) is given by 

w 

T = S. C (£. ) D 2 F. (86) 

v j=A p jv v j 

where D 2v is the transformation matrix from 2-system to v-system. The 
total magnetic actuator force acting on the payload assembly (in v-system) 
is 

F v = D 2v £ F (87) 

j=A J 

The expression for T v in equation (86) is used in the VA rotational 
equation (44) . For the "coarse" assembly, the total reactive force F r 
(in 2"-system) acting on the lateral gimbal is 


F = -C„ f 
r F o 


+ F - e 2 F 
ca z 


rd 


( 88 ) 


where f Q denotes the magnetic actuator output force vector corresponding 
to the command force vector f. The reaction torque on the lateral gimbal 
is given by 


T = T 
r ca 


C„ f 
T o 


+ C (r 
p cm 2 


+ V 


(89) 


The torque acting on the complete "coarse" assembly is given by 


C 


M n - d 02 I CD 02 Pi + D 12 r d + r cm2 + r^) F r 


+ T 

ca 



+ T 


RCJ 


+ T. 


( 90 ) 


in the 0-system, where T dpt and T, denote reaction control jet torque 

KLJ Q 

and disturbance torque. 


30 



These expressions for F^, T^, and are used in the "coarse" 

assembly equations of motion, equations (19), (29), and (36). 

SENSOR MODELS 

Attitude of the shuttle is controlled by vernier reaction control jets. 
In the simulation, however, only one thruster firing is included in order 
to simulate the "worst case" reaction control jet transient. Therefore, 
the shuttle attitude measurement system is not necessary for the simula- 
tion. "Coarse" gimbal positions are measured by sine-cosine resolvers. 

Axial and radial centering errors at the corresponding magnetic actuator 
stations and at the roll motor station are measured by noncontacting 
proximity sensors. Outputs of the proximity sensors are assumed to be 
contaminated with Gaussian zero-mean white noise. The payload attitude 
and rate measurement system consists of three rate gyros and sun sensors/ 
star trackers. The gyros are assumed to be low noise, dry- tuned types, 
and sun sensors are assumed to be high-resolution, monochromatic light 
type. Three rate and attitude error measurements (one per axis) are 
assumed to be available for feeding into the fine-pointing control system. 
Figure 4 shows a single axis block diagram of the payload attitude and 
rate measurement system. Each gyro output is contaminated with a zero- 
mean colored noise and a Wiener process bias. The sun sensor outputs, 
which are discrete-time, are contaminated with zero-mean white noise 
which represents electronics noise, thermal (shot) noise, and quantiza- 
tion error. A bias may also be present in the sun sensor. Gyro drift 
rate, sun sensor bias, and sun sensor boresight misalignment are important 
sources of low-frequency error. Gyro drift rate can be either compensated 
for or estimated by a state estimator which gives an asymptotically unbiased 
rate estimate. Sun sensor bias errors are usually measureable and repeat- 
able, and can be compensated to a large extent by the sun sensor electronics. 
Any uncompensated sun sensor bias or misalignment will affect the mean 
pointing error (pointing "accuracy") , and not the RMS pointing error 
(pointing "stability"). The accuracy requirements are usually far less 
stringent than the stability requirements. Therefore, small sun sensor 
bias errors can be tolerated, and are assumed to be absent. Single-speed 
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wound rotor resolvers are included in the "coarse" gimbal assembly for 
position readout and control of the gimbal angles. A roll axis resolver 
is used to measure payload roll angle relative to the lateral gimbal. 

PAYLOAD ATTITUDE STATE ESTIMATORS 

Pointing performance of the ASPS, or any other pointing system, depends 
on the payload attitude and rate signals used in the synthesis of the 
control torque. If no estimator is used, it can be shown (ref. 4) that 
the pointing performance of a completely . isolated assembly is severely 
limited.. In fact, the pointing error in this case increases directly with 
controller bandwidth, reciprocal of damping ratio, and measurement noise 
level. Reference 8 contains a detailed treatment of the estimator design 
problem. Two types of estimators have been considered in reference 7, 
which give attitude estimates and asymptotically unbiased rate estimates 
despite gyro. bias. The first type does not utilize the input (control) 
torque information, while the second type of estimator utilizes the input 
torque information in its prediction model. In the case of ASPS, the 
input torque can be accurately measured. However, performance of the 
second type of estimator (Filter 3) is critically dependent on the input . 
noise level, which was not accurately known at the time of writing this 
report. Therefore, the first type of estimator (Filter 1A) , one for 
each axis, has. been included in the present simulation. These estimators 
generate estimates <j> v , 0 v , ip , 1^, 0 , ^ for use in the fine 

pointing, control loop. 

SIMULATION DESCRIPTION 

The mathematical model of ASPS/Space Shuttle system consists of a set 
of nonlinear ordinary differential equations. The nonlinearities include 
all trigonometric nonlinearities, magnetic actuator nonlinearities, gimbal 
and gimbal torquer nonlinearities, etc. The control systems consist of a 
magnetic actuator control system, "coarse" gimbal control system, and a 
space shuttle attitude control system. The magnetic actuator control 
system blends the discrete-time signal A av [eq. (83)] generated 
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using payload attitude and rate estimates, with the continuous signal X^ & 
[eq. (84)] to synthesize the magnetic actuator command- force vector. (The 
signal is generated every T^^ seconds and kept constant for 

seconds, where T^^ is the estimator update interval.) Although this 
control scheme has been assumed in the present simulation, other schemes 
such as multirate sampling (different sampling rates for attitude and 
rate feedback) may also be suitable for ASPS application. However, results 
of investigation of these schemes were not available at the time of writing 
this report. Therefore, the single-rate sampling scheme mentioned above 
has been used in the simulation. 

The "coarse" gimbal control system is continuous time, except whert 
used in the backup gimbal pointing mode. In that case,- the control scheme 
would be similar to the fine-pointing control scheme discussed above. 

The space shuttle attitude is controlled by vernier reaction control 
jets. However, for computing RMS pointing errors, a low bandwidth control 
system using an annular momentum control device (AMCD) or control -moment - 
gyros (CMG's) has been assumed for the space shuttle. This assumption is 
not restrictive (ref. 4) because, in the case of ASPS, the payload assembly 
is almost .completely isolated from the shuttle. However, for computing 
peak pointing errors for the "worst case" reaction control jet transient, 
a control jet force pulse of -222.4 N in the Zo direction is applied, 
producing a disturbance torque of -2487.9 N-m about the Xo axis. The 
duration of the pulse is 0.52 seconds, and occurs when the shuttle angular 
velocity about the Xo axis is 0.00435 degree/second. 

It was found that the ac resistance in the magnetic actuator model 
could be ignored without significant effects, reducing the order of each 
(axial and radial) actuator from 6 to 4. Various options are available 
in the computer program. These options include: inclusion or exclusion 

of payload attitude state estimators in the loop; use of magnetic 
actuators or latched vernier assembly backup mode; inclusion or exclusion 
of attitude measurement noise; proximity sensor noise; worst case shuttle 
vernier jet transient; roll motor radial disturbance force; choice of 
payload following or fixed gimbal modes; etc. Dimension of the state 
vector (excluding state estimators) is 65 for normal operation and 30 
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for the backup pointing mode. In order to keep the number of state variables 
low, gimbal angular rates and centering error rates are assumed to be avail- 
able for control synthesis, instead of generating them using second-order 
networks. 

Table 1 lists the noise sources acting on the system, and their standard 
deviations. All noises are assumed to be white and Gaussian. All the noise 
processes, except the sun sensor noise, are continuous time. For digital 
simulation, however, they must be represented by discrete pulses. This can 
be approximately done by choosing an integration interval T^ <_ T^^, and 
generating random numbers with standard deviations given by continuous 
standard deviations divided by /FT. The ranc ^ om numbers generated are held 
constant in each interval T 4 . This sequence of random amplitude pulses 
has approximately the same power as the continuous white noise. In the 
present simulation, T^ , and T^j were both 0.01 sec. The simulation 
proceeds in the following manner: the differential equations are solved 

for T. seconds until time increment T r ., is reached (the ratio T^..,/T. 
is an integer greater than or equal to 1); new state estimates and vernier 
torque commands are generated and held constant for the next T^.^ . seconds 
(magnetic actuator force commands are continuously generated using vernier 
torque commands and continuous proximity sensor outputs); the differential, 
equations are solved again for T^ seconds, and so forth. All noise 
processes are generated every T^ seconds, except the sun sensor noise, 
which is generated every T^^ seconds. 

Mean and RMS errors are computed by solving the differential equations 
until the statistics have evolved, and then computing means and standard 
deviations of the pointing errors as the solution continues. Sensitivities 
to various noise sources are computed by making all except the corresponding 
standard deviation zero. Two options were included in the computer program 
for solving the differential equations --Euler integration and fifth-order 
Adams method. The latter was found to perform satisfactorily for all 
cases because of its ability to solve stiff differential equations. The 
former may be useful for simpler cases (such as backup "gimbal pointing 
systems") . 


34 



NUMERICAL RESULTS AND DISCUSSION 


In order to evaluate the performance of ASPS, a 3.7-m long, 66-cm 
diameter cylindrical payload having a mass of 270 kg was selected. This 
payload is well within the limits of the present magnetic actuators and 
is one of the two payloads used in reference 3. Table 2 gives the 
parameters of the space shuttle orbiter. "Coarse" gimbal torquer and 
friction parameters are given in tables 3 and 4, and cable parameters 
are given in table 5. Payload assembly parameters are given in table 6. 
Table 7 describes the control system modes and parameters used. Table 1 
gives a list of noise sources and their nominal standard deviations. 

In order to evaluate ASPS performance, it is necessary to include the 

errors made in measurement of VA c.m. location and inertia matrix, since 

the measured quantities are used in control systems design. An error of 

-1 percent was assumed in the measurement of r , and of the VA inertia 
r cmv 

matrix I in all the computer runs. The target pointing angles are 
assumed to be <f> = 45°, 0 = 45°, ij> = 0°, which form one of the worst 
combinations (ref. 4) . The basic computer runs for evaluationg ASPS 
performance are described below: 

1. Response to initial condition offsets. 

2. Peak pointing errors for the worst case vernier reaction control 
jet (VRCJ) disturbance. 

3. RMS error computation. 

4. Sensitivities of pointing errors to various noise sources. 

These computer runs are repeated with the following changes: 

a. Use of backup gimbal pointing system (with payload assembly 
latched to lateral gimbal). 

b. No estimators used. 

c. Cable used for data transmission. 

d. "Coarse" gimbals not in the payload following mode. 

That is, the gimbal servos hold them in their nominal positions. 
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Response to an initial condition offset of 20 arcseconds in VA attitude 
and 2 mm (for each axis) in VA translation is shown in figure 5. Although a 
2.5-arcsecond undershoot occurs, the transient dies down in about 1.5 sec. 

The magnetic actuator centering errors decay much more slowly due to low 
suspension bandwidth. 

The worst-case VRCJ transient was next simulated, as described earlier. 
Table 8 shows peak pointing errors under different conditions. Peak pointing 
error was 0.00875 arcsecond (in in the absence of noise, and 0.01006 

arcsecond when noise was present. When the cable was used for data trans- 
mission, the peak error increased to 0.0236 arcsecond. Figure 6 shows 
typical pointing error and centering errors for this transient when the 
cable was used. All peak centering errors were below 4 mm for all these 
cases. These results were obtained for a magnetic actuator current-loop 
bandwidth of 300 Hz. It was found that variation of current loop bandwidth 
had significant effect on peak pointing errors. Figure 7 shows the variation 
of peak pointing errors with current loop bandwidth. 

Theoretically, peak pointing errors should be zero in the absence of 
noise and errors in measuring 1^ and r cmv - However, in reality, the 
peak pointing error is never zero because of (1) a change in VA roll position 
due to the transient, causing a change in the. location of the vernier 
assembly c.m. with respect to magnetic actuators, (2) magnetic actuator 
nonlinearity, that is, error between commanded and actual forces, (3) finite 
bandwidth of magnetic actuator current loop, and (4) sampling of commanded 
magnetic actuator torques. 

The results described above were obtained using the "payload following" 
gimbal mode. Fixed gimbal mode Was next simulated in which the gimbals 
are held by . the control system in a fixed position relative, to the shuttle. 

In this case, it was found necessary to increase the bandwidth of the VA 
translation control system to 0.13 rad/sec in order to limit the magnetic 
actuator centering error excursions . This resulted in somewhat larger 
(about 15 percent higher) peak pointing errors as shown in table 8. This 
demonstrated the advantage of .the payload following mode. 

Peak pointing errors during the backup gimbal pointing mode (with VA 
latched to lateral gimbal) with noise and estimators present, are also 
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given in table 8. The largest peak error for this mode is 1.618 arcseconds 
in 0^-axis, which represents almost a factor of 400 deterioration compared 
to the "normal operation" case. Figure 8 shows the response of the backup 
system to the worst-case VRCJ transient. 

Root-mean-squared (RMS) pointing errors were next obtained as described 
in the preceding section, for various cases. Table 9 shows the results for 
normal operation and for the backup mode. The most significant source of 
error is proximity sensor noise. Contribution of the attitude measurement 
system noise is very small (for the nominal standard deviations). In order 
to investigate the effect of larger attitude measurement noise, the sun 
sensor noise standard deviation was varied from its nominal value, up to 
2000 times the nominal value, and the RMS pointing errors were obtained. 
Figure 9 shows the resulting variation of RMS pointing errors. When the 
cable was used, the RMS pointing error increased only slightly (table 9) 
although the peak error for the worst case VRCJ transient almost 
doubled. 

In order to investigate the effect of roll motor anomaly torque, a 
sinusoidal radial force with a peak of 0.0622 N and a frequency of 200 Hz 
(twice the carrier frequency) was applied at the roll motor station W. 

As discussed earlier, because of the high frequency, a 5-sec time response 
was computed for investigating this effect. All other noise sources were 
held at zero. The resulting RMS pointing errors were (1.12, 0.231, 0.275) x 
10' 4 arcsecond for <p v , 0^, and \/j , respectively, which are insignificant 
compared to the effect -of proximity sensor noise. 

CONCLUSIONS 

A detailed nonlinear mathematical model was developed for the ASPS/ 
Space Shuttle system. All trigonometric nonlinearities as well as nonlineari 
ties in various components such as magnetic actuators, "coarse" gimbals, and 
torquers were incorporated. Control laws and payload attitude state esti- 
mators were designed for solar/stellar pointing missions. A computer 
program was developed to simulate the system, and the pointing performance 
was evaluated for various conditions for normal operation and for the backup 


37 



"gimbal pointing" mode. The RMS errors for normdl operation were (0.9949, 
1.007, 2.844) x 10“ 3 arcsecond for the three axes, and were 0.05359 and 
0.05731 arcsecond in pitch and yaw axes for the backup pointing mode. The 
most significant source of error for the noise parameters considered was 
the proximity sensor noise. Peak pointing errors for the worst case vernier 
reaction control jet transient were 0.00301, 0.00398, and 0.0101 arcsecond 
for normal operation and 1.485 and 1.618 arcseconds for the backup mode. 

Use of a cable for data transmission under normal operation was found to 
cause about a 100 percent increase in peak pointing errors, but minimal 
increases in RMS errors. It can be concluded on the basis of these results 
that ASPS can provide very high quality pointing performance in the presence 
of sensor and actuator noise and disturbances. 
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Measurements Used Standard Deviation 
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Table 2. Parameters of Space Shuttle Orbiter. 


Mass, kg 86176.2 

Inertias, kg-m 2 


I 

XX 

7,882,725 

I 

yy 

1,037,200 

T 

"zz 

8,082,030 

I 

xy 

6,779 

I 

xz 

2,711 

I 

yz 

-178,968 


Crew station location in (X , Y q , Z q ) system, m...(0, -16.002, 0.594) 
Location of point P g in (X q , Y q , Z q ) system, m...(0, 2.2, 0.75) 



Table 3. "Coarse" gimbal parameters. 


Elevation Gimbal 


Lateral Gimbal 


Mass, kg 31.1 103.5 

Inertias (x,y,z), kg-m 2 0.575, 0.575, 0.228 3.84, 3.84, 6.5 


Center of lateral 

gimbal axis in (0, 0, 0.21) 

1' -system (r^) , m 

Location of lateral 
gimbal center 0 in 

2'-system, relative (0, 0, 0.027) 

to lateral gimbal 
c.m. (r £ ) , m 


Location of lateral 

gimbal c.m. in. 2'- (0, 0, 0.208) 

system (r ) , m 
J cm2 


Location of elevation 

gimbal c.m. in 1'- (0, 0, 0.1) 

system (r ) , m 
J cm]/ 
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Table 4. Gimbal torquer and friction parameters. 


Rest stiffness (a^) , N-m/rad 

Peak or Coulomb friction 
torque (T^ ) , N-m 
o 

Percent ripple (p ) 

Ripple frequency (f ) , Hz 
Cogging frequency (f ) , Hz 
Peak cogging torque (T C0 g 0 ) , N-m 


6.77 

0.0316 

2 

48 

96 

0.16 



Table 5. Cable parameters. 


Cable translation spring constants nc . n ^ 

(radial, axial), N/m ’ 

Cable bias forces (radial, axial), N 0.01405, 0 


Cable torsional spring constant (Z)., 
N-m/rad 


0.005 


00445 


Cable bias torque (Z) , N-m 


0.001 



Table 6. Payload (vernier) assembly parameters. 


Mass , kg 

270 

Inertias (x,y,z), kg-m 2 

319, 319, 15 

Location of c.m. in v-system 
(x,y,z), m 

0.07, 0.07, 1.83 


Radius of PMP, m 


0.365 



Nominal Location 

Damping Natural Frequency. 10 of Root due to 

Type of Control ratio, p Nominal Range Integrator 
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Table 8. Peak pointing error. 


Normal 

Operation 


Backup Mode 
(VA latched) 


Peak pointing error, arcsecond 



<f> 

T v 

0 

V 

ip 

r v 

No noise present 

0.00256 

0.00363 

0.00875 

Noise and estimator 
present 

0.00301 

0.00398 

0.0101 

Noise, estimator, and 
cable present 

0.00758 

0.00707 

0.0236 

Noise and estimator 
present, "payload 
following" mode 
not used 

0.00334 

0.00453 

0.0116 

Noise and estimator 
present 

1.485 

1.618 
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Pointing error, arcsecond 
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Figure 3. Definition of position vectors. 




Figure 4. Payload attitude measurement system. 










Figure 5. Response to initial condition offset 
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Figure 6. Response to worst case VRCJ disturbance (with cable) 
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